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Ahhotaiihh. We study the dependence of volatility on the stock price in 
the stochastic volatility framework on the example of the Heston model. To 
be more specific, we consider the conditional expectation of variance (square 
of volatility) under fixed stock price return as a function of the return and 
time. The behavior of this function depends on the initial stock price return 
distribution density. In particular, we show that the graph of the conditional 
expectation of variance is convex downwards near the mean value of the stock 
price return. For the Gaussian distribution this effect is strong, but it weakens 
and becomes negligible as the decay of distribution at infinity slows down. 



1. Introduction 

Stochastic volatility (SV) models are quite popular in recent decades due to a 
need for reliable quantitative analysis of market data. The most popular ones are 
the Heston pi], Stein-Stein [20], Schoble-Zhu [18], Hull- White [32] and Scott [19] 
models. We refer for reviews to [T3], [H], [7]. The main reason for introducing the 
SV models is to find a realistic alternative approach to option pricing to capture the 
time varying nature of the volatility assumed to be constant in the Black-Scholes 
approach. 

Nevertheless, SV models can be used for investigation of another properties of 
financial markets. For example, in [4] the time-dependent probability distribution 
of stock price returns was studied. While returns are readily known from a financial 
data, variance (square of the stock-price volatility) is not given directly so it acts as 
a hidden stochastic variable. In [4] the joint probability density function of returns 
and variance was found, then the integration over variance was performed and the 
probability distribution function of returns unconditional on variance was obtained. 
The latter PDF can be directly compared with the Dow- Jones data for the 20- 
years period of 1982 - 2001 and an excellent agreement was found. The tails of 
the PDF decay slower than the log-normal distribution predicts (the so-called "fat- 
tails" effect). 

Technically our paper is connected with [4]. However, we study the dependence of 
the variance on fixed returns, thus, we estimate hidden stochastic variable through 
the variable that can be easily obtained from financial data. The result strongly 
depends on initial distribution of returns and variance. It is natural that the 
distributions change their shape with time. In particular, we show that for Gaussian 
initial distribution of returns the expectation of variance demonstrates the convexity 
downwards near the mean value of returns. 
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2. General formulas for the conditional expectation and variance 

Let us consider the stochastic differential equation system: 

m dF t = Adt + adW Y , dV t = Bdt + XdW 2 , 

( ' F = f, V =v, t > 0, / e K, v e K, 

where W(t) = (Wi(t), is a two-dimensional standard Wiener process, A — 

A(t,F t ,V t ), B = B(t,F t ,V t ), a = a(t,F t ,V t ), A = \(t,F t ,V t ) are prescribed 
functions. 

The joint probability density P(t,f,v) of random values Ft and V t obeys the 
Fokker-Plank equation fe.g..|17|) 

BP B B 1 B 2 1 B 2 

<2 > m—h w - 1 < bp » + \h { - 2p) + \w {>?p) 

with initial condition 

(3) P(0J,v)=P (/,»), 

determined by initial distributions of Pj and V t . 

If P(t,f,v) is known, one can find P(T4|Pt = /), which is the conditional 
expectation of value Vt at a fixed F t at the moment i. This value can be found 
by the following formula (see, [3]): 

fr t r\ y P(t, f, v )dv 

(4) E(v t \Ft = f)= um J - £ - L) . 

L^ + oo J { LL) P{tJ,v)dv 

Let us also define the variance of Vt at a fixed F t as 

f, r T ,v 2 P(t,f,v)dv 

(5) Var (V t |P t = /)= lim T P u f W, ~ ^ (W = /) • 

In this both formulae the improper integrals from numerator and denominator 
are assumed to converge. The assumption imposes a restriction on the coefficients 
A, B, a, A. 

Note that if we choose Pq(/,«) = 5(v — vo(f))g(f), where vo(f) and g(f) are 
arbitrary smooth functions, then E (Vt\F t = f) |^=o = Vo(f)- 

For some classes of systems ([1]) the conditional expectation V(t, f) was found in 
|17|.[T]. [2] within an absolutely different context. 

Let us remark that sometimes it is easier to find the Fourier transform of P(t, /, v) 
function over /, v variables, than the function itself. We will get formula allowing 
to express E (Vt \F t = f) in terms of Fourier transform of P(t, /, v) and will apply it 
for finding an average variance of the stock price, which depends on known return 
rate. 

Proposition 1. Let P(i,/i,£) be the Fourier transform of function P(t,f,v) over 
(f,v) variables, which is the solution of problem @ ; ([3]), and both integrals from Q 
converge. Assume that P(t,/i,0) and d^P(t,fi, 0) are decreasing over fi at infinity 
faster than any power. Then E (Vt\F t — f) and Var (Vt\F t = f) determined by Q 
and ([5]) can be found as 

6 E(Vt\F t =f)= !L^4-±1ZLJ!}lLLL (>0,/eK, 

' F^[P(t,n,0)](t,f) 
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(7) 

Var(vlF (F^M^^^O^^F-H^^O^F^t^^O)] 

whereF^ 1 andF^ 1 mean the inverse Fourier transforms over fi and£, respectively. 
The proof is a simple exercise in the Fourier analysis. 

3. Example: the Heston model 

Of course, there is no explicit formula for the joint probability density function 
P(t, f, v) for arbitrary system |T]). We will consider a particular, but important case 
of the Heston model 

(8) df t = (a-^dt + ^TtdWu 

(9) dv t = -7(«t - 6)dt + k^/v~ t dW 2 . 
Here a, 7, k, 9 are arbitrary positive constants. 

Equation ((9| describes the process that in financial literature is called Cox- 
Ingersoll-Ross (CIR) process, and in mathematical statistics — the Feller process 
[7] , [6] . In [6] it is shown that this equation has a nonnegative solution for t G [0, +00) 
when 276* > k 2 . 

The first equation describes a return f t on the stock price, in assumption that 
the stock price itself obeys a geometric Brownian motion with stochastic volatility. 
The second equation describes the square of volatility a\ = Vt ■ 

The Fokker-Planck equation (J2J for the joint density function P(t, /, v) of return 
ft and variance vt takes here the following form: 

dP(t,f,v) nu . , , , , a . l2 .dP(t,f,v) , f v \ 8P(t,f,v) , 



7 P(t, /, v) + ( 7 (« - 0) + fc 2 ) ^|P + (I - a) 



(10) 



dt ' v ' ■" ' y ' v ; ' <9v V 2 / 9/ 

k 2 vd 2 P(t,f,v) vd 2 P(t,f,v) 



2 <% 2 2 df 2 

Now we can choose different initial distributions for return and variance. Note that 
it is natural to assume that initially the variance does not depend on return. 
Below we denote E (v t \ft — f) as V(t, f) for short. 

The function P(t, the Fourier transform of P over (f,v), satisfies the 

equation 

8P %^ + \ 0* + ^ + 2 7 C + *k 2 e) dP{t d ^ ] + i (j8 + ^a) P (t, M , = 0. 
The first-order PDE ([TT]) can be integrated, the solution has the following form: 

(12) p, t ^ _ c t^± p (' (»7f + »)->■■ - ggj j 

2iarctan( - fc ^ +t7 ) 



(13) *F /x,-t + 



v 7 ? 
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where q = —ik 2 [i + k 2 ji 2 + 7 2 , F is an arbitrary differentiable function of two 
variables. 

3.1. The uniform initial distribution of returns. We begin with the simplest 
and almost trivial case. Let us assume that initially the rate of return is distributed 
uniformly in the interval (— L, L), (L = const > 0), and volatility is equal to some 
constant a > 0. Then the initial joint density distribution of f t and v t is 

(14) P{0,f,v) = ±-6{v-a). 

To simplify further calculations we will exclude randomness for t = 0, i.e. we will 
assume a = 0. 

The respective initial condition for the Fourier transform is 

(15) P(p,p,Z) = ~6(n). 
The solution of problem (fTTj) , (fT5l) takes the form 

■yet 

(16) P(t, t x,Z) = -S( f i) 1 1 



L~ yr ~' \ (2 7 eT* + ik 2 £ ( e f* - l)f ) 
It is easy to calculate that 

(17) P(t lM) 0) = p(/i), 9 € P(t, M ,0) = p(M)^(e-^-l). 



Finally from © and we get 

(18) E(v t \f t = f)=e(l 

Ok 2 9 

(19) Var(v t \f t = f) = — (l - e^f . 

It is evident that here there is no dependence on / and the result is the same as 
we could obtain from calculation of mathematical expectation and variance of v t 
from equation ©. 

3.2. The Gaussian initial distribution of returns. Let us assume that initially 
rate of return is distributed according to the Gaussian law. Then we have the 
following initial condition: 

(20) P(0, /, v) = ^ e - m2f2 S(v), m>0. 

When a = 0, the Fourier transform of initial data over (/, v) is P (0, fj,, £) = e f- 1 . 

Solution of the problem (fTTj) , (|20|) takes the form: 
(21) 

-k 2 t + i"f \\\ ~^ 



cosh I — y k 2 fi(fi — i) + 7 2 — i arctan 



yjk 2 \i{\i — i) + 7 2 
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We see that P(t, /x, £) exponentially decreases over /x. That is why we can use 
formula ([6]) and obtain (after cumbersome transformations) the following integral 
expression: 

(22) V(tJ)=2 7 6 r 



where 



*(t, ft, f) = e ^ 



(A cosh (f) + 2 7 sinh(f )) 
sinh (^) 



(Acosh(f)+2 7 sinh(f))' 
A 2 = fc 2 (V + 1) +4 7 2 . 

Let us remark that if a =/= 0, we can also get a similar formula, but it will be 
more cumbersome. 

The limit case asm->oo for (12011 is 



(23) P(Q,f,v) = 6(f)5(v), m>0. 

For this case the formula (|22l) modifies as follows: the exponential factor in the 
expression for \1/ takes the form e lfJ, ^~ ta \ 

3.3. "Fat-tails" initial distribution of returns. Integral formula, analogous to 
(|22p can be obtained for initial distributions intermediate between uniform and 
Gaussian ones. For example, as initial distribution we can take 

P(0, /, v) = K(l + m 2 f) q 8(v), m > 0, q < 0, 

with an appropriate constant K . Exact formula for the Fourier transform P(t, ^i, £) 
can be found for q = — |, — n, n £ N. For all these cases P(t, fi, £,) decays as \fi\ — > oo 
sufficiently fast and Proposition [T] can be applied for calculation of V(t, f). 



For example, for q = — 1 the difference with (|21l) is only in the multiplier e i- 2 : 
it should be changed to 

with the Heaviside function H. 

3.4. Convexity downward of the volatility curve and asymptotic behavior 
for small time. It turns out that if in the Heston model the average volatility is 
considered as a function of the rate on return, we will observe a deflection of the 
plot. The effect appears in numerical calculation of both integrals in (f2"2"| with 
the use of standard algorithms. The numerical calculation of the integrals over an 
infinite interval is based on the QUADPACK routine QAGI [TB], where the entire 
infinite integration range is first transformed to the segment [0,1]. For example, 
Fig. 1 presents the graph of function V(t, f) at three consequent moments of time 
for the following values of parameters: 7=l,fc = l,^ = l,Q! = l ) m = l. 

This behavior of the volatility plot can be studied by analytical methods as 
well. Indeed, let us fix rate of return /. Then from (f22|) by expansion of integrand 
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functions into formal series as t — > up to the forth component and by further 
term-wise integration (series converge at least for small / and m) we will get that 
(24) 



v(t, /) = jet - \-i 2 et 2 + ± 7 ( 7 2 + 2/ 2 TO 4 fc 2 - /TO 2 fc 2 - m 2 fc 2 ) t 3 

-±70 (87 fc 2 / 2 TO 4 - 4 (7 + 4rn 2 a) k 2 m 2 f - 4 (7 + a) m 2 k 2 + 7 3 ) t 4 + 0{t 5 ). 

Let us justify a possibility to expand V(t, /) into the Taylor series. We should prove 
that both integrals in the numerator and denominator of (|22[) can be differentiated 
with respect to t. Indeed, let t £ fl t — (— r, r), < t < 00. It can be readily shown 
that both integrands in (|22l) . $ and ty, are continuous with respect to /i and t on 
Ixfij, the derivatives of any order <9" $, <9" "J, n — 0,1,... are also continuous 
on M x f2 t . Moreover, |9™<I'| can be estimated from above by c\ ■ e~ C2AI , 

with positive constants c\ and C2. Therefore J R <9™<I>d/x and L d^^d^i converge 
uniformly on VL t . Thus, according to the classical theorem of calculus the numerator 
and denominator in (|22[) can be differentiated on tt t under the integral sign. Since 
for 

J Q(t,n,f)dfi\ t = = J d?Q(t,fx,f)\ t = d», 

Q = or $, the Taylor coefficients in the expansion of J R Q(t, /x, /) d/i can be 
obtained by integration of the respective coefficient of Q(t, /) of the Taylor series 
in t with respect to /i. The latter integrals can be explicitly calculated. This gives 
expansion (|24"]) . 

Hence for t —> we find that 

V(t,f) ~ ±70t 3 (l-7t)m 4 fc 2 / 2 -(I am 2 t+\(l - 7 1)) i 3 7 9m 2 k 2 f+ 



i (7 1 - (l 3 - at)) 9 jt 3 m 2 k 2 + 7 6 1 - | 7 2 (9 i 2 + 1 7 3 ? - ± 



is a quadratic trinomial over / with a minimum in point / = 4 4^ ra 2( t 1 _T ) f t ) 1 1 f° r t > 0. 

The effect holds for initial "fat-tails" power initial distributions as well. 
Nevertheless, this effect weakens as the decay of the distribution at infinity becomes 
slower. 

Fig. 2 presents the function V(t, f) for three consequent moments of time for the 
initial distribution of return p(f) given by formula 

P(f) = ~ 1 i 1 1 2 ' 

7T 1 + / 

The values of parameters are 7 = 10, Ar = 1, 6 = 0.1, a = 10. It seems that the 
curves are strait lines, but the analysis of numerical values shows that the deflection 
still persists near the mean value of return. Acting as in the case of the Gaussian 
initial distribution one can find the Taylor expansion of V(t, f) as t — > 0, 

v(t, f) = yet - W - .e^^e + ^ MfM t . + 0{t% 

where we denote by a polynomial of order k with respect to /. We do not 
write down these polynomial, let us only note that ^j^fjj^ ~ (js — an< ^ 

~ (2 7 2 - as 00. 
It is very interesting to study the asymptotic behaviour of V(t,f) as |/| — > 00 
and i — >■ 00. We do not dwell here on this quite delicate question at all and reserve 
it for future research. Some hints can be found in [4], [9], [10] . 
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3.5. Modifications of the Heston model. Let us analyze the situation when the 
coefficient 7 from equation ([9]) depends on time. For some interesting cases of this 
dependence one can find the Fourier transform of P(t, /, v) and formula for V(t, /). 
For example, if we set 7 = jct, then we get a Brownian bridge- like equation (see, 
|15| describing square of volatility behavior with start at vq — a > and end at 
Vt = b > 0. Here the solution will be represented in terms of integrals of Bessel 
functions and the solution is cumbersome. 

It may seem that the described approach, which helps to find the conditional 
expectation of volatility under fixed returns in the Heston model, can be successfully 
applied in other variations of this model. This is true when initial rate on return 
has a uniform distribution. However, this situation is trivial, because the answer 
does not contain / and is equal to the expectation of return obtained from the 
second equation of model. In the case of non-uniform initial distribution of return 
(for instance, Gaussian) formula ([6]) may be non-applicable, even when explicit 
expression for P(t, n, £) can be found. The cause is that P(t, n, £) increases as 
\fi\ — > 00. For example, if we replace equation @ with 

(25) dv t = -7(« t - 6)dt + kdW 2 , 7, 0, k > 0, 

under initial data (|2U)) . a = 0, we will get 

~ . ^ ( k 2 t 4 k 2 t , fet k 2 t 1 \ 2 [9 \ \ 

whence it follows that the coefficient of /i 4 in exponent power is positive when t is 
positive. This means that integrals from ([6]) are divergent. 

4. Possible application 

Basing on our results one can introduce a rule for estimation of the company's 
rating based on stock prices. The natural presumption is that company's rating 
increases when return on assets increases and volatility decreases. Hence for 
estimation of the company's rating one can use (very rough) index R(t, f) — 
f/V(t,f), where V(t,f) is calculated by formula (|22l) . Figs. 3 and 4 shows the 
plot function R(t, f) for three consequent time points for Gaussian and power 
distributions, respectively. Parameters as in Figs. 1 and 2. We can see that in 
the Gaussian case the index does not rise monotonically with return. 
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5. Conclusion and further work 

In this article we obtain an estimate of volatility given rate on return data in 
the frame of the Heston model. This problem has been solved by calculation of 
the average volatility under a fixed rate on return and under the supplementary 
condition on initial distribution of return and volatility. Namely, different cases 
of initial distribution of returns have been studied: uniform, Gaussian and "fat- 
tails" distributions, intermediate between them. We revealed that the graph of tthe 
averaged volatility is convex downwards near the mean value of the stock price 
return for the Gaussian initial distribution and for certain distributions decreasing 
at infinity slower then the Gaussian one (for which we succeed to find the Fourier 
transform of the joint probability density of return and variance explicitly) . For the 
Gaussian distribution this effect is strong, but it weakens and becomes negligible 
as the decay of distribution at infinity slows down. 

Let us note that our formulas can be obtained in a different way, using the 
well-known expression for the joint characteristic function of the log-return and 
the variance in the Heston model [TT| (in the correlated case) . This expression was 
obtained exploiting the linearity of the coefficients in the respective PDE, in other 
word, the fact that the Heston model is affine [5]. Nevertheless, this way is not 
convenient for our purpose, since it requires an additional integration. 

Formulas for the conditional variance at fixed return V(t, /) are obtained in 
the present work in the integral form, we compute the integrals numerically using 
standard algorithms and study asymptotics of the formulas for small time. The 
questions on analysis of the formulas for larger t and / and on the asymptotics of 
V(t, /) as |/| — > oo and t — > oo are open. Moreover, the dependence of the averaged 
variance on the properties of the initial distribution of returns has to be studied in 
general case, not only for separate examples, as it was done here. 
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